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ABSTRACT 

We present an analysis of the near-infrared continuum emission from the cir- 
cumstellar gas disks of Be stars using a radiative transfer code for a parametrized 
version of the viscous decretion disk modeL This isothermal gas model creates 
predicted images that we use to estimate the HWHM emission radius along the 
major axis of the projected disk and the spatially integrated flux excess at wave- 
lengths of 1.7, 2.1, 4.8, 9, and 18 /im. We discuss in detail the effect of the disk 
base density, inclination angle, stellar effective temperature, and other physical 
parameters on the derived disk sizes and color excesses. We calculate color excess 
estimates relative to the stellar V^-band flux for a sample of 130 Be stars using 
photometry from 2MASS and the AKARI infrared camera all-sky survey. The 
color excess relations from our models make a good match of the observed color 
excesses of Be stars. We also present our results on the projected size of the disk 
as a function of wavelength for the classical Be star ( Tauri, and we show that 
the model predictions are consistent with interferometric observations in the H, 
K', and 12 /im bands. 

Subject headings: stars: emission-line. Be — circumstellar matter — infrared: 
stars — stars: individual {( Tau) 



1. Introduction 



Be stars are rapidly rotating B-type stars that at times host a circumstellar gas disk 
( jPorter fc Riviniusll2003l ) . They exhibit various observational attributes such as photospheric 
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non-radial pulsations, hydrogen and iron emission lines that are formed in the disk, infrared, 
millimeter, and radio continuum excess emission from the disk, intrinsic linear polarization 
from electron scatter ing in the disk, and, in an increasing; number of cases , the presence of a 
compact companion (jStruvelll93ll : ISlettebaklll988l : iPorter fc Riviniusll2003[ ). The circumstel- 
lar ionized gas around Be stars gives rise to a large infrared e xcess emission that is seen in 
their spectral energy distribution. Observations conducted by IWoolf et al.l (Il970l ) at 10 /xm 
were the first to suggest that the IR excess of Be stars is due to free- free and bound -free emis- 
sion and not due to thermal dust emission. An early study by iGehrz et al.l ( 1l974j ) confirmed 
that the large IR emission detected from 33 Be stars originates in the free-free emission 
only and that thermal dust emission fails to reproduce the observations at 20 /im. The IR 



excess from ionized dis 



fICote fc Waters 



1987 



i gas was subsequently investigated yith s pace-b ased data frorn . IRA S 



Waters et all 119871 : iDoughertv et al.lll994l ) ISO flWaters et al.ll2000h 



and MSX (IClarke et al.ll2005l ). and with ground-based mid-IR (IRinehart et al.l Il999l ) and 
near-IR data from the 2MASS survey (jZhang et al.ll2005l ). The near- and mid-IR excess flux 
emission is known to increase with wavelength and to dominate over the stellar flux at long 
wavelengths. Observations probing the circumstellar disk structure at many wavelengths 
and across all Be spectral types are necessary to understand the complex physical processes 
involving disk formation. 

Models of the disk continuum emission have become increasingly comprehensive in recent 
years. Some o f the first models that assu r ned a simple geometry and gas density distrib ution 
JWatersI Il986l : iKastner fc Mazzalil Il989l : lOoughertv et al.l 1 1994 iRinehart et aD Il999[ ) were 
successful in fitting the color excesses in many cases. A significant irnprovenient came with 
the d evelopment of the viscous decretion disk model (ILee et al.lll99ll : |Porteiill999l : lOkazaki 
200ll ) in which the gas orbits with nearly-Keplerian velocity, but with a small radial outflow 
motion caused by turbulent viscosity. These models characterize the disk density with a 
radial power law, p oc R~^'^, and a vertical extension set by hydrostatic equilibrium. The 
IR-excesses predicted are similar to those observed, but often require a density exponent with 
smaller absol ute value and /or a non-isothermal temperature distribution in order to match 
observations (jPorterl Il999l . |2003| ). Recent models include radiation processes in the disk 



and determine fully con sistent temperature and density distributions (ICarciofi fc Bjorkman 



20061 : I Jones et al.l l2008l ) . and these are remarkably successful in explaining the continuum 



and line emission in some Be stars (ICarciofi et al.ll2006l : iHalonen et al.ll2008l ). 



The advent of optical long baseline inte rferometry has now led to the direct resolu- 



tion of Be star disks in both emission lines (|Ouirrenbach et al.l 119971 : IStee fc Bittai 



2007 



Meilland et al 



2007 



Tvcner et allbood: Meilland et al.|[2007l : Carciofi et al.ll2009 ) and in the IR- continuum (j 



2001 



Gies et al. 



200 



^ . These observations offer us the means to explore the geom- 



etry of Be disks and to measure the disk gas temperature. Since the IR-excess flux depends 
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mainly on the projected size of the disk and the gas surface brightness (from the tempera- 
ture dependent source function), a comparison of the angular size and flux excess provides 
a way to investigate the disk-to-star temperature ratio and (in principle) to study spatial 
temperature variations in the disk. 

Our goal in this paper is to explore the IR-excess flux predictions for a parametrized 
version of the isothermal, viscous disk r aodel that we h ave developed to study CHARA 
Array i^'-band observations of Be disks (IGies et al.l 120071 ). We calculate the near-IR color 
excesses for various representative stellar and disk paramete rs and compare these with new 
IR measurements from the AKARI satellite (llta et al.ll2010l ). In §2, we describe the model 
in detail and discuss how the different physical parameters of the model influence the color 
excess. In §3, we discuss the color-color diagrams and compare our results to the AKARI 
all-sky survey observations of Be stars. Finally in §4, we consider the example of the Be 
star ( Tau, and we show that the predicted variation of circumstellar disk angular size with 
wavelength is consistent with the available multi-wavelength observations. 



Be Star Circumstellar Disk Models 



Our original model (iGies et al.ll2007l ) created images of Be stars and their disks for com- 
parison with i^'-band interferometry from the GSU CHARA Array. In order to investigate 
the disk emission over a range of wavelengths in the near and mid-IR, we have extended the 
model to calculate the disk flux distribution at 1.66 (//-band), 2.13 (i^'-band), 4.8, 9, and 
18 fim. Our code is a realization of an isothermal viscous decretion disk model. In brief, the 
gas density distribution is given by 



p{R, Z) = poR-"" exp 



Z 



2 V^(^) 



where R and Z are the radial and the vertical cylindrical coordinates, respectively, in units 
of stellar radii, po is the gas base density, and n is the radial density exponent. H{R) is the 
disk vertical scale height defined by 



H{R) 



Vk 



(2) 



where Cg is the sound speed (in turn dependent on the disk gas temperature) and Vk is the 
Keplerian velocity at the stellar equator. Our multi-wavelength disk model has four phys- 
ical parameters: the base density po, the density exponent n, the disk-to-star temperature 
ratio Td/Tcs, and the outer boundary disk radius -Rout- There are also two observational 
parameters: the wavelength A and the inclination angle i. 
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The code follows the method of iHummel fc VranckenI (120001 ) and solves the equation of 
transfer along a grid of sight-li nes around the s tar and disk. The result is a spatial image 
of the system (see examples in lGies et al.l 120071 ). We assume source functions equal to the 
Planck functions for the temperatures of the star and disk. The disk optical depth in the 
near-IR is mainly due to free-free and bound-free processes, and it can be expressed using 
an incremental step ds along a projected rectilinear coordinate grid as follows 



dr = C{\,Td)p{R, Zf ds 



(3) 



where the coefficient C(A,Td) is given by equation (5) in iDougherty et al.l (119941 ). This 
coefficient includes terms for the Gaunt factors for bound-free and free-free emission. We 



set th ese Gaunt factors by interpolating in wavelength in the tables from IWaters fc Lamers 
(Il984j ). For simplicity, we estimate the ionizat ion equilibrium arid Gau nt factors throughout 
the d isk for two idealized cases (adopted from lLamers &: Waterslll984j and lWaters &: Lamers 
19841 ): (1) a hot plasma with ionized H, singly ionized He, and doubly ionized C, N, O for 
disk temperatures above 15000 K, and (2) a warm plasma consisting of ionized H, neutral 
He, and singly ionized C, N, and O for lower temperatures. Note that we ignore the electron 
scattering flux in the model, since this flux source is proportional to the stellar flux and 
hence is relatively small at longer wavelengths. 

Here we focus on two results of the calculation, the flux excess and apparent disk size. 
We determine from the derived image the total monochromatic flux of the star and disk, 
Ftotai, and then we estimate the net disk contribution as Fd = Ftotai — Fg, i.e., the net flux 
relative to the unobscured star. We show this quantity in a magnitude form of flux excess 
as E*(y* — m\) = 2.51og(l + F^/Fs). This notation is based upon the assumption that the 
disk contribution is negligible in the optical \^-band, so that we can write the flux excess 
as how many magnitudes brighter the system appears at longer wavelengths where the disk 
shines brightly. This latter assumption should be taken with caution. In fact, there is ample 
evidence from the V^-band brightening of Be stars during their active phases that the disks do 
contribute to the V^-band (perhaps by as much as 50%) in those cases with dense and large 
circumstellar disks. Consequently, by referring the colors to a V magnitude that is brighter 
than that for the star alone, the observed color excesses may be lower than the calculated 
color excesses in some cases. The star superscript is used here to differentiate this type of 
"reddening" o r color-excess from the kind normally associated with interstellar extinction. 
Waters! (119861 ) expresses the monochromatic flux excess as Z = -Ftotai/-^s; so in our notation 
E*(y* — nix) = 2.5 log Z. We also use the calculated spatial image of the star plus disk to 
find the HWHM of the emission envelope along the projected major axis of the disk, and we 
define an effective, observational disk radius Rd/Rs as the ratio of the angular HWHM to 
the angular stellar radius. 
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We begin by showing our results on the flux excess and disk size for a default model, 
and then we show how changes in the physical and observational parameters affect the 
results. In the default model, we assume that the central star is an early-type star with 
effective temperature Tes = 30 kK, radius Rs/Rq = 10, and mass Ms/Mq = 15.5. The 
infrared excesses derived from observations suggest that the power- law density exponent of 



the c ircumstellar disk falls in the range n ~ 2.0 to 3.5 (ICote fc Waters! 119871 : IWaters et al. 



987f). so ^ye assu med n = 3 here, a value consistent with our prior interferometric results 



(iGies et al.l 120071 ). The other adopted parameters for the default model are a disk-to-star 
temperature ratio Td/Tes = 2/3, an outer boundary disk radius Rout/Rs = 21.4, and an 
inclination ang le i = 45° . Th ese parameters are selected from an earlier model for the 



Be star 7 Gas (iGies et al.l 120071 ). and the outer boundary, for example, corresponds to the 
Roche radius of this binary system. Our results for the color excess and disk radius are 
listed in Table 1 as a function of waveband and base density for this model and several 
others described below. 

We show in Figure 1 how the color excesses E*(y* — K), E*{V* — Q/im), and E*(y* — 
IS/^m) vary as a function of the disk base density po- The disk flux excess is highly dependent 
on density po and wavelength A. The color excesses at low densities are insignificant because 
the disk is optically thin in the continuum and the stellar photospheric flux dominates. 
As the density increases, the color excess resulting from the disk emission becomes more 
important and dominates at longer wavelengths. In fact, the excess emission at 18 /im is 
higher than at 9 pm and at 2.13 /im because the optically thick-thin boundary of the disk 
becomes bigger and the excess flux larger at longer wavelengths. 

The relationship between the 18 /im color excess and apparent disk emission radius for 
the default model is shown as a solid line in Figure 2 (for T^/Tcg = 2/3). We also show 
a number of other models where different parameters are varied in turn. The case with a 
higher inclination angle i = 80° is plotted with a dashed line. As the inclination increases, the 
projected disk surface area decreases and hence the disk flux excess also declines. However, 
at higher inclination, a ray through the outer part of the disk encounters significant density 
over a longer path, and the increased optical depth causes the effective radius to appear 
larger. Gonsequently, as the inclination increases, points on the default model curve are 
shifted to lower color excess and higher effective radius. 

Ghanges in the other parameters have less influence on the size - color excess relation 
shown in Figure 2. For example, for the case of a cooler Be star with an effective temperature 
of Teff = 15 kK (shown by the dashed-dotted line), the resulting color excess and disk radius 
are almost identical to that for the default case (Teff = 30 kK; solid line) at all values of the 
disk base density. This is due to the temperature dependence of the optical depth coefficient, 
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C{\^Ti) oc . The radius of the optically thick-thin boundary will vary with this optical 
depth term, and the source function is proportional to temperature in the Rayleigh- Jeans 
part of the spectrum. Thus, the disk emission flux will vary as the product of projected area 
and source function, or ~ i?^ S" ~ ^d^^'^Y Td, which is approximately constant, all other 
parameters being equal. 



If we adopt a smaller value of the outer radius of the disk [Rout/Rs = 14.6; shown by the 
dash-triple dotted line), we see that there is no difference between this case and the default 
case {Rout/Rs = 21.6) at low densities. It is only at very high density that the truncation of 
the outer disk leads to a slight decline in the color excess. 

Finally, the dotted line shows the influence of our choice of disk temperature. In the 
default mo del, we adopt an isoth e rmal disk with a temperature of = 2/3 T^s following the 
example of Hummel fc VranckenI ( 2000 ). We compare in Figure 2 the color excess for two disk 
temperatures, 1/2 T^s and 2/3 Tgg. The source function varies approximately linearly with 
disk temperature Td, so a drop of 25% in disk temperature will create a comparable decrease 
in emission flux (for a given disk radius). This decrease of ~ 0.3 mag in E*{V* — ISfim) is 
seen in Figure 2 at the high density end where the disk flux dominates. This suggest that 
coordinated interferometric and IR excess observations are potentially an important means 
to study the disk temperature properties, especially for Be stars with dense and large disks. 



Our results show many similarities to the color excesses derived by iDougherty et al. 



(I1994J ) from a much simpler model. They conflned the emitting gas to a wedge-like disk 
with a density law dependent only on the dist ance to the star, and the y solved the radiative 



transfer problem for just i = and i = 90°. 'Dougherty et al.l (119941 ) also investigated the 



dependence of the near-IR color excess on the model parameters, but they present relative 
color excesses between adjacent near-IR bands {CE[J — K) and CE{K — L); see their Fig. 7) 
rather than referencing the excess to the stellar \^-band flux as we do here. Nevertheless, the 
agreement is reasonable, and, for example, their expression for the disk optically thick-thin 
boundary as a power law function of wavelength a nd base density (t heir e g. 7) is similar 
to our results for Rd/Rs in Table 1. Furthermore, IDougherty et al.l (119941 ) also presented 
instructive results for a range in the density power law exponent, with values greater and 
less than the n = 3 value assumed here (they used the symbol (3 for the density exponent). 
Their work showed that a smaller n yields a more spatially extended disk and consequently 
larger color excesses. 

Recent models for Be star disks present detaile d calcu lations of the gas temperature 
as a function of disk position. ICarciofl &: BjorkmanI (120061 ) solve for the temperature dis- 
tribution in the disk through a Monte-Carlo treatment of radiative transfer, and they find 
that gas temperature reaches Td ~ 0.6 T^s in optically thin parts of the disk in one rep- 
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resentative model. ISigut et al.l ( 120091 ) present results for a grid of models that maintain 
hydrostatic equilibrium with the spatial variations in temperature, and they find that the 
density averaged temperature is ~ O.GTeff, however, this average temperature declines 
with increasing disk density (perhaps by 20 — 30%; see their Fig. 8). Their models indicate 
that the temperatures tend to be cooler in the denser regions closer to the star and in the 
mid-plane. The IR-ex cesses calculated frorn such models are similar to those for isothermal 
models. For example, ICarciofi fc BjorkmanI (120061 ) show that the predicted flux excesses for 
different inclination angles are almost the same for the isothermal and non-isothermal model 
they study (see the upper panel of their Fig. 10). However, since the disk gas temperatures 
decline with increasing density, we suspect that a plot like Figure 2 of color excess and radius 
for such fully consistent, non-isothermal models would look similar to our — 2/ ST^eff curve 
at small excess (low density), but would tend towards the = l/2Tcs curve at large excess 
(high density). 



AKARI IR Fluxes of Be Stars 



With the recent release of the AKARI/IRC mid- infrared all-sky survey ( llshihara et al. 



2010l ). we now have the opportunity to compare the observed and model flux excesses at 9 
and 18 fim. We took as our sample 130 Be stars from the work of lDougherty et al.l (Il994[ ) that 
have reliable estimates of interstel lar redden i ng Ej B—V). We then collected V, Kg, m(9 /im), 
and m(18 /zm) magnitudes from llta et al.l ( l2010l ) for each targe t. The V magnitudes were 
taken from the SIMBAD database. Kg magnitudes from 2MASS ( ISkrutskie et al.ll2006l ). and 
m(9/im) and m(18/xm) magnitudes from AKARI. These last two are based on the Vega 
magnitude scale, where a model spectrum b y R. L. Kurucz defines the flux zero point as a 
function of wavelength ( iTanabe et al.ll2008l ). The photomet ry we use here was colle cted at 
different times, and since Be stars are inherently variable (jPorter &: Riviniud |2003[ ). some 
scatter must be expected in the results because the flux excesses will change with disk density 
variations. We also collected stellar effective temperatures for each target from the work of 
Fremat et al.l (|2005[ ) in order to estimate the intrinsic stellar colors. 



Our goal is to determine a flux excess by comparing the observed and intrinsic stellar 
colors of the targets. Using the same magnitude notation given in the previous section, we 
determine the near-IR flux excesses by 



V -mx- E{B - V") X (3.10 - Rx) - {V - mA)(Kurucz) 



wher e the ratio of interstellar extinction to reddening is R\ = A\/E{B — V) ( iFitzpatrick 
19991 ) and (V — mA) (Kurucz) is the intrinsic stellar color derived from monochromatic sam- 
pling of flux ratios of model spectra with the Vega spectrum from R. L. Kurucz. These 
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model spectra are also from Kurucz atmospheres for solar metallicity, gravity \ogg = 4.0, 
and a microturbulent velocity of 2 km s~^ (parameters appropriate for main sequence B- 
stars). The intrinsic colors are listed in Table 2 as a function of effective temperature Tgfj, 
and we found that they made a reliable match to the AKARI colors o f B-stars with known 
Teff from interferometry and bolometric luminosity (jCode et al.l 119761 ) . The derived color 
excesses for the Be stars are listed in Table 3. The typical errors in the AKARI magnitudes 
are ±0.05 mag, but they can be larger for the 2MASS Kg magnitudes since many of the 
Be stars are bright and their magnitudes were determined from the the wings of the point 
spread function. We caution that the intrinsic colors from Table 2 may not be appropriate 
in some cases because rotational gravity darkening will make stars with large inclination 
appear redder (although the color difference may be negated if opaque disk gas blocks the 
cooler equatorial zones from view). We have ignored these complications because they are 
difficult to estimate accurately and because these colors are not too sensitive to temperature 
for hot stars. 

We plot the results in two color-color diagrams in Figures 3 and 4. These show the flux 
excesses at 9 and 18 /xm, respectively, as a function of the i^^-band excess. Also shown in 
these figures are plots of our model near-IR excesses (§2) for the default model and cases with 
differing stellar temperature and disk inclination. All these models make similar predictions 
about color-color excesses, and they appear to match the observations well, especially if 
allowance is made for the slight negative shift in observed colors caused by our neglect of disk 
fiux in the l^-band. There are three objects plotted in Figures 3 and 4 with E*{V* — K) ^ 1.0 
that fall well below the predicted tre nds. All three ar e binary stars in which t he companion 
is a K -giant (HD 45910 = AX Mon. lEhas et al.l ll997l: HD 50123 = HZ CMa. Isterken et al. 



19941 : HD 50820, iGinestet fc I2OO2I ). and we suspect that their relatively large 



brightness in Kg is due to the fiux from the cool giant companion. Otherwise, the overall 
agreement suggests that the viscous decretion disk model provides a satisfactory description 
of the near-IR fiux excesses. 

Inspection of Figures 3 and 4 shows that at the high density limit, the color-color excess 
diagrams assume a linear form. This part of the relation occurs when the disk fiux dominates 
over the stellar fiux and the colors become those of the disk. Thus, we expect that all Be 
stars will appear with approximately the same near-IR color for those cases with sufficiently 
dense disks. As an example, we show in Figure 5 a near-IR color-magnitude diagram for 
all those sampl e stars with Hippar cos parallax data yielding absolute magnitude errors less 
than 0.5 mag (Ivan LeeuwenI 120071 ). This plot shows the interstellar extinction corrected, 
absolute Kg magnitude versus an interstellar reddening corrected, color index Kg — m(9 fim) 
(shown as plus signs). Each of these is connected by a dotted line to a color and magnitude 
coordinate corresponding to one with the derived color excess removed (nominally for the 
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star itself; shown as diamonds) . Also plotted is the zero-age main- sequence (left dashed line) 
from the work of iLejeune fc Schaererl (120011 ) that was formed from the theoretical {T^g, V) 
track and the colors in Table 2 for T^s = 10 — 30 kK. From Figure 3, the asymptotic form of 
the excess relation is 

E*{V* -9 i2m)^ E^iV - Ks) + 1.35 

and we also plot the main-sequence translated in color by this expression and brighter in Kg 
by 1 magnitude (see Fig. 1) to represent the approximate positions of the dense disk case 
(right dashed line). We see that the Be stars appear over a range in color between the no 
disk and strong disk cases, with many of the Be stars having colors close to the strong disk 
lir nit. This diagra r n is si milar in appearance to that for the Be stars discovered in the LMC 
by lBonanos et al.l (120101 ). who present a plot in the (J — [3.6yLim], J) plane. 



4. Angular Size of the Disk of C, Tau 



Meilland et al.l (120091 ) measured some of the first Be disk diameters at 8 and 12 /xm 
using the VLTI and MIDI instrument. Their results suggested that the disk sizes do not 
increase with wavelength as predicted by simple models. For example, they found that the 
angular size of the disk of the Be star a Ara was approximately constant between 2 and 
12 yum, and they speculated that the disks may be truncated by the tidal effects of a binary 
companion. Our models and those of [Dougherty et al.l (Il994j ) suggest that the color excesses 
and effective radii could reach finite limits in those cases with high disk density and a small 
outer boundary. 

Here we demonstrate that at least for one case, C, Tau, the angular size variation with 
wavelength is consistent with model predictions. The Be star Tau (HD 37202, HR 1910 , 
HIP 26451) is a frequently observed target with a strong IR-excess (ITouhami et al.ll2010l ). 
The Ha emission line in its spectrum shows cyclic V/R variations on a timescale of few 
years , which are explained by the presence o f a one-armed oscillation in its circumstellar 
disk (jOkazaki et al.l 120021 : ICarciofi et al.l 120091). The star is the primary in a spectroscopic 
binary with an orbital period of P = 133 d (IRuzdjak et al.ll2009l). The system is com posed 



of a IIMq primary Be type star and a 1.3M0 seco ndarv star (jFloauet et a.1 



eral interferometric studies o f the Ha emission line (iQuirrenbach et al 



11997 



19891). Sev- 



Tycner et al. 



20041 ) and the IR-continuum Jcies et al.ll2007l : ICarciofi et allbood : iMeilland et al.ll2009l ) have 



resolved th e circumstellar disk a round the primary Be star. A recent CHARA Array inves- 
tigation by lSchaefer et al.l (120101 ) shows that the disk of C Tau is viewed almost edge-on and 
that the disk may be precessing with the V/R cycle. 



We used our model to estimate Rd/Rs for C Tau over the wavelength range of 1.7 
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18 yum. We adopted the stellar parameters from iGies et al.l (120071 ): mass Ms = II.2M0, 
radius Rs = 5.5Rq, effective temperature Teg = 19 kK, and parallax vr = 7.82 mas. The 
disk parameters assumed are po = 1.4 x 10^^° g cm~^, n = 2.9, i = 87°, and T^/Tcs = 2/3 
(based on fits of recent K' observations with the CHARA Array). The outer disk boundary 
was set at the binary Roche radius, i?out = -RRoche = 146i?0. We then determined the disk 
effective radius over the wavelength grid, and our results for Rd/Rs are plotted in Figure 6. 

In order to test our model predictions, we collected recent interferometric measurements 
of the angular size of ( Tau. We start with the two CHARA Array r esults. The weighted 

p[) is 1.61 ± 0.05 mas. 



average of the if-band, Gaussian FWHM from ISchaefer et al.l (120 



and the corresponding i^^'-band value is 1.79 ± 0.07 mas (IGies et al 



2007f). T he target was 



also observed at longer wavelengths with VLTI/MIDI by lMeilland et al.l (120091 ). They found 
an upper limit of FWHM less than 2.6 mas at 8 fim, but they resolved the disk at 12 fim 
and found a FWHM = 5.7 ± 2.2 mas. These estimates are over-plotted for comparison on 
the theoretical dotted-line shown in Figure 6 by assuming a stellar, limb darkened, angular 
diameter of Old = 0.40 ± 0.04 mas. We find that the size of the emitting region does 
increase with increasing wavelength in a manner mostly consistent with the observations. 
The 8 fim upper limit falls slightly below the model curve in Figure 6, but we suspect that 
this difference is marginal given the fact that the VLTI/MIDI observations were not made 
at an optimal position angle for the disk's orientation in the sky. Thus, we suggest that the 
viscous decretion disk predictions about disk size as a function of wavelength pass the test 
for the case of C Tau. However, it is certainly possible that disk truncation effects will be 
more important in binary Be stars with smaller semimajor axes. 

We are currently completing a survey of bright, northern sky. Be stars with the CHARA 
Array in the K'-hand, and we plan to use our model to fit both the angular diameters and 
flux excesses. Since the flux excess depends on the angular size and source function, we will 
be able to infer the source function and hence disk temperature for a diverse sample of Be 
stars. We also intend to comp are the H and . fi^^-ba nd diameters for a subset of Be stars with 
CHARA MIRC observations JSchaefer et allboioh . 
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Table 1. IR Color Excesses for Different Viscous Disk Models 













PO (g cm 










Parameter 


1.0 X 10 


12 


5.1 


X 10-12 1.0 X 10"" 


3.1 X 10- 


11 


8.1 X 10" 


in 

1.0 X 10" " 


in 

2.0 X 10" " 






n 


= 3, 


i = 45°, Tefi 


= 30 kK, Rout 


= 21.4i?s 


, = 


= 2/3 Tefi 






E*(V -H) .... 


0.02 






0.03 


0.06 


0.29 




0.91 


1.08 


1.65 


E''{V* -K) .... 


0.02 






0.05 


0.11 


0.59 




1.38 


1.57 


2.19 


E*{V* -4.8 ixm) 


0.02 






0.11 


0.31 


1.04 




1.86 


2.06 


2.69 


E*lv -9 /im) . 


0.04 






0.43 


0.91 


1.86 




2.76 


2.97 


3.61 


E*{V* - 18 Atm) 


0.13 






1.00 


1.61 


2.59 




3.51 


3.71 


4.33 


RrtlR^iH) 


1.23 






1.24 


1.25 


1.41 




2.10 


2.28 


2.97 


Rri/RAK] 


1.23 






1.24 


1.26 


1.58 




2.39 


2.61 


3.40 




1.24 






1.27 


1.39 


2.17 




3.15 


3.44 


4.56 


Rj/R„(Q um) 


1.24 






1.38 


1.79 


2.76 




4.07 


4.49 


6.05 


B,,/BJ_1H /ini) 


1.26 






1.95 


2.55 


3.92 




5,95 


6.57 


8.82 






11 






= aO kK, Rom 


= 21.1Ks 


, I 'j = 








E*{V -H) .... 


0.02 






0.03 


0.04 


0.11 




0.39 


0.49 


0.91 


E*{V - K) .... 


0.02 






0.03 


0.04 


0.25 




0.72 


0.87 


1.41 


E*{V* - 4.8 ^Jin) 


0.02 






0.05 


0.13 


0.48 




1.11 


1.28 


1.88 


E* {V* - 9 ^lIa) . 


0.03 






0.19 


0.43 


1.11 




1.95 


2.15 


2.78 


E*{V* - 18 Atm) 


0.07 






0.49 


0.89 


1.78 




2.68 


2.88 


3.48 




1.24 






1.26 


1.31 


1.83 




2.77 


2.95 


3.99 




1.24 






1.27 


1.39 


2.14 




3.13 


3.49 


4.56 




1.24 






1.39 


1.80 


2.87 




4.29 


4.62 


6.21 


R^/R,(9 um) 


1.26 






1.77 


2.43 


3.64 




5.61 


6.11 


8.17 


R^/R«(18 am) 


1.37 






2.57 


3.43 


5.33 




8.03 


8.86 


11.97 






n 


= 3 


i = 45", Teff 


= 15 kK,Rout 


= 21.4i?a 


Td = 


= 2/3 Teff 






E*{V -H) .... 


0.02 






0.04 


0.08 


0.43 




1.11 


1.28 


1.86 


E*{V* - K) .... 


0.02 






0.06 


0.15 


0.73 




1.53 


1.73 


2.35 


-4.8 /im) 


0.04 






0.17 


0.46 


1.26 




2.11 


2.31 


2.94 


E*{V* -9 (Um) . 


0.05 






0.49 


1.01 


1.96 




2.86 


3.07 


3.71 


E*{y* - 18 um) 


0.14 






1.04 


1.64 


2.63 




3.54 


3.79 


4.36 


Ra/Rs{H) 


1.23 






1.24 


1.27 


1.64 




2.46 


2.68 


3.51 


Rd/RsiK) 


1.23 






1.25 


1.29 


1.77 




2.62 


2.85 


3.76 


i?d/i?s(4-8 Atm) . 


1.24 






1.30 


1.55 


2.44 




3.59 


3.94 


5.23 


Rd/Rs{9iJ.m) ... 


1.24 






1.52 


2.04 


3.12 




4.63 


5.12 


6.86 


Rd/Rs{lS IJ.m) .. 


1.27 






1.99 


2.62 


4.03 




6.07 


6.71 


9.03 






n 


= 3, 


i = 45", Tefi 


= 30 kK, Rout 


= 21.4i?s 


,Td-- 


= 1/2 Teff 






E*(y* - H) .... 


0.02 






0.03 


0.06 


0.30 




0.91 


1.08 


1.65 


E*\v* -K) .... 


0.02 






0.04 


0.09 


0.45 




1.16 


1.34 


1.94 


E*{V* - 4.8 ^lm) 


0.03 






0.11 


0.31 


1.04 




1.86 


2.06 


2.69 


E* {V* - 9 ^lm) . 


0.04 






0.31 


0.72 


1.61 




2.48 


2.68 


3.33 


E*{V* - 18 ^lm) 


0.09 






0.79 


1.35 


2.32 




3.22 


3.42 


4.05 


Rd/Rs{H) 


1.23 






1.24 


1.25 


1.41 




2.10 


2.28 


2.97 


Rd/RsiK) 


1.23 






1.24 


1.26 


1.55 




2.34 


2.54 


3.31 


Ra/Rs(-i.S Atm) . 


1.24 






1.27 


1.39 


2.17 




3.15 


3.44 


4.56 


Rd/Rs(9tim) ... 


1.24 






1.39 


1.79 


2.76 




4.07 


4.49 


6.06 


i?d//ls(18 /im) .. 


1.26 






1.86 


2.46 


3.75 




5.67 


6.25 


8.48 






n 


= 3, 


j = 45", Tefi 


= 30 kK, Rout 


= 14.6i?s 


, Td -- 


= 2/3 Teff 
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Table 1 — Continued 



PO (g cm ^) 



Parameter 


1.0 X 10-12 


5.1 X 10-12 


1.0 X 10-" 


3.1 X 10-" 


8.1 X 10-" 


1.0 X 10-1° 


2.0 X 10-1° 




E*{V* - H) .... 


0.02 


0.03 


0.06 


0.29 


0.91 


1.08 


1.65 


E*(V* - K) .... 


0.02 


0.05 


0.12 


0.60 


1.38 


1.57 


2.18 


E*{V* -4.8 Mm) 


0.02 


0.11 


0.31 


1.04 


1.86 


2.06 


2.67 


E*lv -9 /iiii) . 


0.05 


0.43 


0.92 


1.86 


2.75 


2.96 


3.57 


E*{y* - 18 tan) 


0.13 


1.00 


1.59 


2.58 


3.47 


3.66 


4.23 


Rd/Rs{H) 


1.23 


1.24 


1.25 


1.42 


2.10 


2.29 


2.97 


Rd/Rs{K) 


1.23 


1.24 


1.27 


1.58 


2.39 


2.60 


3.40 


Rd/Rs{i.8 lira) . 


1.24 


1.27 


1.39 


2.17 


3.15 


3.44 


4.56 


Rd/Rs{^ tiui) ... 


1.24 


1.38 


1.79 


2.76 


4.07 


4.48 


6.05 


Rd/Rs{18nm) .. 


1.26 


1.95 


2.55 


3.92 


5.95 


6.56 


8.82 



Table 2. Adopted Main Sequence Colors 



Teff V -K V -9 urn V -18 urn 
(kK) (mag) (mag) (mag) 



10 


-0.06 


-0.07 


-0.07 


12 


-0.23 


-0.29 


-0.29 


14 


-0.34 


-0.43 


-0.44 


16 


-0.43 


-0.54 


-0.55 


18 


-0.51 


-0.64 


-0.65 


20 


-0.58 


-0.74 


-0.75 


22 


-0.64 


-0.82 


-0.83 


24 


-0.70 


-0.90 


-0.91 


26 


-0.74 


-0.96 


-0.98 


28 


-0.79 


-1.02 


-1.03 


30 


-0.83 


-1.07 


-1.08 



- 16 - 



Table 3. Be Star Color Excesses 



E*{V* - K) E*{y* - 9 //m) E*{V* - 18 ixm) 



nu 


(mag) 


(mag) 


(mag) 


144 


0.05 


0.12 




4180 


0.13 


1.54 


2.44 


5394 


0.54 


1.83 


2.51 


6811 


-0.34 


-0.02 


0.41 


10144 


0.03 


0.29 


0.43 


10516 


0.46 


1.68 


2.68 


18552 


0.11 


0.74 


1.62 


20336 


0.03 


1.47 


2.16 


22192 


0.33 


1.75 


2.62 


23016 


0.03 


0.06 




23302 


-0.07 


0.72 


3.24 


23480 


-0.06 


0.29 


1.72 


23552 


-0.07 


0.53 


1.54 


23630 


0.29 


0.28 


0.97 


23862 


0.10 


0.45 


1.25 


25940 


0.08 


1.16 


1.99 


28497 


0.30 


1.87 


2.51 


29866 


0.03 


0.72 


1.49 


30076 


0.57 


1.63 


2.60 


32343 


0.39 


1.21 


2.43 


32990 


-0.05 


0.07 




32991 


0.46 


1.56 


2.65 


35439 


-0.23 


1.54 


2.12 


36576 


0.70 


1.88 


2.95 


37202 


0.72 


1.71 


2.56 


37490 


-0.15 


0.96 


1.68 


37795 


0.10 


0.51 


1.30 


41335 


0.52 


1.90 




44458 


-0.19 


1.49 


2.21 


45542 


0.00 


0.32 


0.89 


45910 


1.02 


1.79 


2.15 


46860 


0.06 


0.13 




50013 


0.23 


1.49 


2.21 


50123 


0.90 


1.60 


2.08 


50820 


1.01 


0.80 


1.33 


54309 


0.39 


1.00 


1.90 


56014 


-0.04 


0.88 


1.40 


56139 


0.01 


1.28 


2.36 


57150 


0.40 


1.67 


2.52 


57219 


-0.06 


0.12 




58155 


-0.04 


0.00 




58343 


-0.25 




2.19 


58715 


-0.03 


0.47 


1.17 


60606 


0.52 


1.41 


2.57 


60855 


0.75 




3.09 
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Table 3 — Continued 



E*{V*-K) E*{V* -9 urn) E* {V* - 18 iJ.m) 



HD 


(mag) 


(mag) 


(mag) 


63462 


0.28 




1.94 


65875 


0.39 


1.73 


2.64 


66194 


0.64 


1.52 


2.69 


68980 


0.48 


1.62 


2.45 


71510 


0.01 


0.09 




72067 


-0.02 


1.08 


2.03 


75311 


-0.09 


0.19 


1.05 


77320 


0.39 


1.16 


1.41 


79621 


0.03 


0.02 




83953 


0.50 




2.37 


86612 


0.10 


1.57 


2.58 


88661 


0.54 


1.78 


2.69 


91120 


0.06 


0.46 


1.21 


91465 


0.32 


1.51 


2.23 


105435 


0.29 


1.64 


2.49 


107348 


0.03 


0.55 




109387 


0.29 


1.42 


2.30 


110432 


0.26 


1.82 


2.50 


113120 


-0.44 


1.33 


2.09 


120324 


0.07 


1.69 


2.24 


120991 


0.29 


0.07 


1.03 


121847 


-0.03 


0.00 




124367 


0.42 


1.56 


2.50 


127972 


0.10 


1.13 


1.87 


137387 


0.07 


0.45 




138749 


-0.01 


0.07 


0.20 


142184 


0.09 


0.54 


1.18 


142926 


0.11 


0.48 


1.05 


148184 


0.49 


1.93 


2.71 


153261 


0.83 


2.09 


2.92 


156325 


0.00 


0.12 




157042 


0.20 


1.34 


2.02 


158427 


0.64 




2.48 


158643 


0.39 




4.40 


164284 


-0.13 


0.05 


0.59 


164447 


-0.07 


0.12 




167128 


0.03 


0.62 


1.34 


168797 


-0.02 


0.24 




170235 


-0.13 


0.07 




171780 


0.29 


0.12 




173370 


0.04 


0.26 


0.28 


173948 


0.12 


-0.03 


0.23 


174237 


0.21 


0.75 




175869 


0.03 


0.06 




178175 


0.29 




2.40 
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Table 3 — Continued 





E*{V* - K) 


E*{V* - 9 ixm) 


E*{V* - 18 iJtxa) 


HD 


(mag) 


(mag) 


(mag) 


183362 


0.55 


1.77 




185037 


0.04 


0.61 




187567 


0.56 


1.78 


2.63 


187811 


0.33 


0.35 


1.38 


189687 


0.04 


0.33 


1.46 


191610 


0.25 


0.14 


0.78 


192044 


0.12 


0.81 


1.37 


193911 


0.00 


0.51 


1.23 


194244 


0.04 


0.24 




194335 


0.19 


1.34 


2.41 


195554 


0.03 


0.17 




196712 


-0.04 


0.53 




197419 


0.28 


0.88 




198183 


-0.09 


-0.08 


0.28 


199218 


0.08 


0.84 




200120 


0.45 


1.69 


2.13 


200310 


-0.03 


0.26 


1.79 


202904 


0.04 


1.62 


2.50 


203025 


-0.32 


-0.12 




203467 


0.58 


1.91 


3.47 


205551 


-0.14 


0.06 




205637 


0.14 


0.69 


1.37 


208057 


—0.06 


0.03 




208682 


-0.28 


-0.14 




209014 


0.04 


0.62 


1.32 


209409 


0.22 


1.10 


1.94 


209522 


-0.01 


-0.06 




210129 


-0.09 


1.24 


1.99 


212076 


0.41 


1.74 


2.48 


212571 


-0.41 


0.51 


0.91 


214168 


0.25 


0.62 




214748 


-0.02 


0.14 


0.99 


216057 


-0.07 


0.01 




216200 


0.25 


0.63 




217050 


0.28 


1.65 


2.58 


217543 


-0.07 


0.33 




217675 


-0.12 


0.04 


0.44 


217891 


-0.05 


0.93 


1.58 


224544 


0.00 


-0.07 




224559 


0.10 


1.52 


2.15 
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Fig. 1. — Plots of the variation in color excess as function of disk base density for the default 
model parameters {n = 3, i = 45°, Tcs = 30 kK, i?out = 21ARs, Td = 2/3Tee). The excess 
emission increases with increasing disk base density po and wavelength. 
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Fig. 2. — Plots of the ratio of disk HWHM radius to stellar radius as a function of the flux 
excess at 18 //m. The sohd line illustrates the relationship for the default model while the 
other line styles show the results found by changing one of the model parameters (indicated 
in the legend) . 
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-1 I .... I .... I .... I .... I .... ^ 

-0.5 0.0 0.5 1.0 1.5 2.0 

E*{V*-K) (mag) 



Fig. 3. — A color excess - color excess diagram showing the 9 /xm flux excess as a function 
of the Ks-hand excess. The solid line shows our results using the default model input 
parameters, the dashed line represent a model with a high- inclination disk {i = 80°), and the 
dashed-dotted line shows the case of a cooler Be star (Tcs = 15 kK). The plus signs indicate 
the observed excesses of Be stars derived from 2MASS and AKARI/IRC photometry. 
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E*{V*-K) (mag) 



Fig. 4. — A color excess - color excess diagram showing the 18 //m flux excess as a function 
of the Kg-hand excess (in the same format as Fig. 3) . 
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Fig. 5. — A color - magnitude diagram showing the absolute Kg magnitude corrected for 
interstellar extinction as a function of the interstellar reddening corrected, color index — 

m(9 yum) for our Be star sample (shown as plus signs). Each target point is connected by a 
dotted line to a diamond representing the color and magnitude of the B star alone. The left 
dashed line represents the zero-age main-sequence for stars with Tpfj = 10 — 30 kK, and this 
also is shown translated in color and magnitude for a dense disk as the dashed line on the 
right hand side. 
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Fig. 6. — The disk-to-stellar radius ratio as a function of wavelength for our model of the Be 
star ( Tau (dotted line). Also shown are the measured sizes from interferometry (or upper 
limit on the size for the 8 /im observation) based on an adopted stellar angular diameter of 
0.40 mas. 



